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ONSET OF ANOMALOUS DIFFUSION IN ELECTRON -BOMBARDMENT 

ION THRUSTOR 
by Allan J. Cohen 
Lewis Research Center 

SUMMARY 

Noise measurements and variations of ion production energy per beam ion with mag- 
netic field indicate an onset of anomalous diffusion in the electron-bombardment ion 
thrustor. Application of the stability theory of the Kadomtsev and Nedospasov kind 
yielded an expression for the critical magnetic field for the onset of anomalous diffusion 
and an expression for the oscillatory frequency associated with the anomalous diffusion. 
The derived behavior of the critical magnetic field is similar to an experimentally de- 
rived scaling law for the ion thrustor. Calculated values of the critical magnetic field 
for mercury and argon propellants agreed well with those obtained from experimental 
values determined from curves of ion production energy per beam ion. Expressions for 
ion production energy per beam ion are given for both the classical and the anomalous 
regions. A theoretical value of the oscillatory frequency was calculated to be in the same 
frequency regime as an experimental value of a device dominated by a similar type of 
anomalous diffusion. 


INTRODUCTION 

The electron-bombardment ion thrustor was developed for use in electric propulsion. 
This engine has been the subject of many investigations, and its operation and perfor- 
mance are well documented (refs. 1 to 3). The thrustor consists of an ion source in the 
form of a plasma discharge exhausting into an electrostatic accelerator system. In the 
ion source, a plasma is created through the electron bombardment of neutral propellant 
atoms. In the electron bombardment process, electrons are conducted across an axial 
magnetic field by the application of a radial electric field between a cathode and an anode. 
The magnetic field serves the purpose of containing the electrons long enough so that they 
can be utilized most effectively in the electron-bombardment process. 



The process of conduction and diffusion of charged particles in a magnetic field has 
recently been the subject of many investigations (ref. 4). As a result of these studies, 
an anomalously high value of conduction or diffusion (here and afterwards referred to as 
anomalous diffusion) has been found in many of the experiments. This anomalous diffu- 
sion results in the conduction or diffusion of charged particles at rates much higher than 
can be expected from classical calculations. Furthermore, this anomalous diffusion fre- 
quently tends to start beyond some critical value of magnetic field. Included in these in- 
vestigations is a wide range of geometries of electric and magnetic fields as well as a 
wide range of plasma densities. 

There are two indications of anomalous diffusion in the plasma geometry associated 
with an ion thrustor. First, an onset of noise has been found beyond some critical value 
of magnetic field (refs. 5 and 6); this behavior is typical of the onset of anomalous diffu- 
sion. Second, curves of ion production energy per beam ion (here and afterwards re- 
ferred to as eV/ion) for the ion thrustor show an expected decline with an increasing mag- 
netic field up to approximately this same critical value. Beyond this value of the magnetic 
field is a transition into either a leveling off region or an increasing region, a result, as 
it will be shown, that cannot be explained classically. 

In this report an application of the stability calculations first introduced by 
Kadomtsev and Nedospasov (ref. 7) is made to the plasma source of the ion thrustor. 

(The theory was first applied to the positive column of a high density plasma discharge 
and later successfully reworked to a low density regime by Guest and Simon (ref. 8). ) In 
this theory an oscillatory type of perturbation is introduced into the number density and 
electric potential of the plasma source. The condition that the growth of the perturbation 
becomes undamped, and thus continually increases with time, yields a critical magnetic 
field. In the study herein, the theory of Kadomtsev and Nedospasov is modified to apply 
to a type of plasma source having a different relative orientation of electric and magnetic 
fields from that to which it was originally applied. Also, the theory is applied to a lower 
plasma pressure regime. The results of this application yield a critical magnetic field 
and a frequency of oscillation, both associated with the onset of anomalous diffusion. The 
critical magnetic field can be compared with values obtained from eV/ion curves of the 
ion thrustor, and the frequency of oscillation obtained can be compared with experiments 
of a plasma source of similar design. 

The curves of eV/ion against the magnetic field of the ion thrustor behave in a clas- 
sical manner up to a critical magnetic field and thereafter show either a constant value 
or an increase with increasing magnetic field. The classical behavior, as well as the 
anomalous behavior, of the variation of eV/ion with magnetic field is described in this 
report. 
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ANALYSIS 


Anomalous Diffusion Mechanism 

The mechanism of anomalous diffusion introduced in the Kadomtsev and Nedospasov 
instability theory is based on a macroscopic distortion where the perturbed number den- 
sity and potential have a 6 -direction dependence. As a result of this dependence, anoma- 
lous diffusion is caused by a Hall effect which results in increased charged particle cur- 
rents in the direction transverse to the magnetic field. In the plasma column to which 
this theory was first applied, the electric field had the same direction (the z-direction) 
as the longitudinal magnetic field. In the plasma source of the ion thrustor, the applied 
electric field which can be considered as the driving force is in the r-direction perpen- 
dicular to the longitudinal magnetic field. There is an indication that with this configura- 
tion the same type of anomalous diffusion mechanism is applicable. First, an analysis 
in reference 9 considers a similar type mechanism as that of Kadomtsev and Nedospasov 
and applies it to a plasma source with an axial magnetic field and an inwardly directed 
radial electric field. The result of this analysis is the prediction of critical magnetic 
fields which are in agreement with experimental results. Further, an experimental study 
(ref. 10) of a low density plasma source with magnetic and electric fields perpendicular 
to each other in a cylindrical geometry resulted in the finding of anomalous diffusion. In 
this experiment, oscillatory variations as well as a 6 -direction dependence of number 
density were found and shown to cause the anomalous diffusion. In the present report the 
Kadomtsev and Nedospasov type of anomalous diffusion mechanism is applied with the re- 
sult that good agreement is found between derived and experimental critical magnetic 
fields. The mechanism is further applied in this report to the anomalous eV/ion curves 
of the ion thrustor. 


Unperturbed Calculation 


The basic geometry of the plasma source in the ion thrustor is shown in sketch (a). 
The region which is analyzed is a central portion of the plasma source, far from both 
ends (the long cylinder approximation). The normal magnetic field of an ion thrustor is 
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strong enough to affect greatly the motion of the electrons but has negligible effect on the 
ions. For the electrons the condition 


CO T » 1 

e e 

is satisfied. (Symbols are defined in appendix A. ) Ions are produced throughout the re- 
gion of analysis by electron bombardment of neutral atoms. The equation of motion for 
electrons is (ref. 8) 


kT — 

— — Vn = — vXB+ — E - — (1) 

m e n m e m e T e 

The electrons are tied to the magnetic lines of force and diffuse across the chamber as a 
result of neutral-electron collisions. The electric -field force on the ions is directed ra- 
dially inward, but the ions can only reach the cathode by collisions in which there is a 
loss of angular momentum. For the electrons the equations of motion derived from equa- 
tion (1) become 


3n 


n e v r, e 


= -D 


J-,e 


3r 


3V1 

dr I 


+ n e^l, e ~ 1 " W - T 


D i,e ^e M l,e Sv} 
e e \ r 30 e r 30 , 


( 2 ) 


n e v 0, e 


T /p ^ ^ 

r 30 r 30 / + e e l e 3r eMi ’ e 3r, 


(3) 


nv = 0 
e z, e 


(4) 


where the Hall currents are taken into account. The motion of the electrons in the 
z-direction of the thrustor is assumed to be zero. For the ions the equations of motion 
become 


3n. 


n.v i = [-D, . — -- n./x, . — 
i r,i \ l,i 0r i^l,i ^ 


3V 

9r , 


(5) 


n i v 0, i 


_ Pj -> i 8n i _ n M l,i dv} 
< r 30 * r 30 ) 


( 6 ) 


n i v z, i = 0 


(7) 
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The Hall currents for ions, because of the weak magnetic field typically used in ion thrus- 
tors, are negligible. The motion of the ions in the z-direction should be small since the 
analysis is restricted to the central portion of the ion thrustor, and the axial potential 
gradient there is small. Appendix B is included to consider currents in the z-direction. 
The continuity equations can be written for the present analysis as 


9n 

— +v - n e f e = Zn e (8) 

8n. _ 

— 1 + V • n.v. = Zn (9) 

11 e 


where Z is the number of ionization events per electron per unit time. In the system, 
plasma neutrality is assumed so that 


n e = n i = n 


Substituting equations (2) to (7) into equations (8) and (9), with the assumption of plasma 
neutrality, yields the following equations: 


— " D i i 

at » *1 


i -L (r + ^1- is, t |i i- LiX) + i- i. fn») 

r Sr\ 3r/ 2 dQ 2\ ’ x r 0r \ 3r/ 2 90 \ M / 


= Zn (10) 




r 3r \ 3r/ f 2 ^2 ,e r dr \ dr/ r 2 d 0 \ 30/ 


"eVi.e, 

ten dV _ 3n 3V) 


r 1 

(3r 30 30 dr) 

1 " 


= Zn (11) 


Next, the normal unperturbed state (n = ng) can be deduced with regard to the follow- 
ing conditions from equations (10) and (11). The steady state is recovered by taking 
an Q /at as zero. Taking n^ and Vq (the unperturbed potential V) as being independent 
of 0 results in dn^/ 30 and 3Vq/ 30 being equal to zero. The resultant equations are 
then 
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(12) 


D l, j a L 9n 0 ) j 9 


rn. 


r 3r \ 3r 


r 3r 


dVr 


3r 


= Zn, 


0 


and 


P l, e 3 L 8n o\ , \ . i 


rn r 


r 3r \ 3r 


r 3r 


dV r 


3r 


= Zn r 


Now eliminating the second term between equations (12) and (13) yields 

- (D i, i^l, e + D l, e M l, i> ~ ^ ^ 7 ) = Z(Mi > e + i )n 0 


which is of the form 


(13) 


(14) 


3 2 n n « 3n n 
2 + 1 + An n = 0 

3r 2 r 9r 


(15) 


where A is a constant (see appendix B for the constant including z-currents). The solu- 
tion of equation (15) is 


n 0 (r) = n|,J 0 (ar 0 r) 


(16) 


where a^r is the argument of the zero-order Bessel function, and n^ and c*q are con- • 
stants that can be solved for by the insertion of boundary conditions. Taken as a boundary 
condition, np(a) = 0 (ref. 7) yields 


a 


0 = 


2.4 

a 


(17) 


where 2.4 is the first zero of the zero-order Bessel function, and r * a is the radius 
of the anode. 


Stability Calculations 

Critical magnetic field. - To consider the stability of the plasma source of the ion 
thrustor, a small perturbation is added to both the unperturbed number density nQ and 


6 



(18) 


the potential Vq such that 


n = n Q + nj 


and 


v = v 0 + v i (19) 

Substituting these quantities into equations (10) and (11) and ignoring products of perturba- 
tions yield the following equations: 


3nj 

at 


D 




1 _0_ 
r 3r 



M l,i 


i _a_ 

r Sr 


rn. 


^0 

3r 


+ rn, 



= Zn, 


( 20 ) 


3n n 

— - D i e 

3t ’ e 


13 / 9n i\ 1 02n i 

+ M, p 

"1 3 L 3V o 8V A . "o d2y i 

r I + 

r 1 3r Hr J r 2 _ 

_ r 0r V 0r / r 2 30 2 _ 

1 , e 


jit , _ /3n n 3V 1 3n., 3V n \ 

e e% _e l_0 __1 _ __1 _0 a z (21) 

r \dr d6 30 3r J L 


In equations (20) and (21) the unperturbed equations have been subtracted out and the 
quantities 3^/30 and 3 Vq/ 30 are again taken equal to zero. Next, the time dependence 
and angular dependence of the perturbed quantities are taken as 

n x = ni (r)e iwt+ime (22) 


and 


Vj = V 1 (r)e ia,t+im ^ 

When these are substituted into equations (20) and (21), the results are 


(23) 
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i&rn^r) - Zn,(r) - D ± •<- — i 
1 A x,1 lr dr 



3ni(r) 
r» 1 

Jr 3r 

3r 



n 





<l± 

3V g av^rf 

rn, (r) — ^ + rn n 

N 

- ±5? nQV^r) ► 

r dr 

< 

1 3r u 3r 

r J 


> = o (24) 


and 


ton l( r).Zn l( r)-D 1>e ^i 


k 

3n,(r) 
r ± 

[r 3r 

dr 


- -L m 2 ni (r)J 


+ ^1, e^~ T“ 

’ r 3r 


,, 9V o 

rnj(r) + rn f 


3V, (r) 


0r 


3r 


-^ n o v iW 

r z 


CO.T^j. 


e e^x, e 
r 


3n n 3Vn 

— -imV.,(r) - imn, — - 
3r 1 1 3r 


= 0 


(25) 


For the type of instability considered, the fundamental mode m can be taken equal to 
unity. Solving equations (24) and (25) directly is difficult and, instead, an approximate 
technique of Kadomtsev and Nedospasov is used. The technique is to take specific func- 
tions of r that satisfy the boundary conditions and substitute them into the' equations. 
The test functions used are 


n x (r) = n , 1 J 1 (a 1 r) (26) 

and 

Vj(r) = VjJjtejr) (27) 

where Jj(a ^r) is the first-order Bessel function of argument a^r. The boundary condi- 
tions used to evaluate are that n^(a) = 0 and Vj(a) = 0 so that 

<*.= — (28) 

A a 
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where 3. 8 is the first zero of the first-order Bessel function. Substituting equations (26) 
and (27) into equations (24) and (25) and noting that 


n 0 = n 0 J 0 (o! 0 r) 


(29) 


9 Jjtejr) 1 W^r) J^r) _ , 

QfjjjtQ'jr; 


3r 2 


r 0r 


(30) 


A [xJj(x)] = xJ Q (x) 


(31) 


yield 


p 

icon^c^r) - Zn , jJ 1 (o' 1 r) + Dj^ .o; jn , Jj(o' 1 r) 


9V 0 9 V 0 
- Ml, i n W a l r ) — a l - Mi fi Ji(«i r ) — n i 
dr 9r^ 


-i-k 


rJ 0 (“o r> 


0J 1 (a 1 r)‘ 


0r 


J 0 (“o r > , X 


> » 0 (32) 


and 


iwn , 1 J 1 (Q' 1 r) - Zn^J^a^r) + D i > e Q! l n i J l( Q! l r ) 


, 9V o 9 v o , 

+ Mj. ,e n l J 0^1 r ^ — a l + Mi , e J 1 ( Q! 1 r )— “ n l 
9r 0r 2 


- e n b V ’f " ~ 

■ L ’ e u 1 r 0r 


rJ 0 (« 0 r) 


0J 1 (o' 1 r)' 


0r 


J o< a o r > T , x, 

+ J l(“l r > 


ia> e T eMl e ^ T //M ^ T //w ^ , iCt, e T eMi e „t T ^ 9V 0 n 

+ Q! 0 n Q V jJ ^a^r^ j(Q! x r) + —7 n x J ^a- jr) — — = 0 


(33) 


9 



For the purpose of simplification, the following quantities are defined: 


I / “i rJ ° 

Jo 

if- 

Jo 


0V n 

(«t r )Ji(«i r ) — - d(«ir) = L 
1 1 1 0r 1 


a S 2 V 

jrj|(a!jr) - d(ajr) = M 

0r 2 



J 0 (« 0 r ) 


J l(«l r )Md(«l r ) s N 


i C Q' 1 a! 0 J 1 (a 0 r)J 2 (Q' 1 r)d(Q; 1 r) ^ q 
C •'O 


where C is defined as 


J' Q' 1 rJ 2 (o' 1 r)d(a!jr) = C (34) 

Multiplying equations (32) and (33) by Q?^rJ^(a ^r)/C and integrating from 0 to a yield 

» 

[iw - Z + a i - ii lf j(L + M)J — i + (M if jNjVi = 0 (35) 

L n n 


and 


[to - Z + a 2 Dl> e + M e (L + M) + to e T eMi _ e p] 4 + (-M ±J e N + to e r e(li> e Q)Vi = 0 

n n 


(36) 
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The results are two homogeneous algebraic equations. The determinant of these equa- 
tions is set equal to zero. This yields 


- N "i, e [*" - z + “X i ' 1*1, i< L + M)] + 1 1 [to - Z + afD 1( , - ^ t (L + M)] 


N "l,i 


iaj - Z + a e + Mx e (L + M) + i — 
’ ’ B 


= 0 (37) 


where w e T e Mx e = (1/ B ). Taking the real part of equation (37) yields 


± , e + X, i^^I - Z > + “l N W, e D l, i + X, i D l, e } = 0 


(38) 


where co = + ia^ (the real and imaginary parts of co) and w^/B is small with re- 
spect to terms such as „ When o>j goes from negative to positive, the perturbed 

number density and electric potential grow without limit, and hence the condition that 
= 0 is the critical condition. The term ojj is thus set equal to zero, and the Einstein 
diffusion relation, / u = (eD/kT), is used. The equation 


.z(Xl± 

\ kT 


le D l /— + — V 

kT i / ’ ’ \ kT e kT i/ 


(39) 


is then obtained from equation (38). Noting that for this analysis the conditions 


! eD . 

±z 1 » x > e 


kTj 


kT, 


and 


e » e 


kT. 


kT, 


are satisfied gives the result that 


Z = 


“l D JL, e 


(40) 


where a ^ = (3. 8/a). The following substitutions can be made into equation (40): 
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Z - V’iV 


D i . * i v „r? 

J-j G g Gj H L 


V- „ = (or v A ) 
e, n n ' e, n e x 


r L = 


m v, . 
e i, e 

eB 


where ( > implies averaging over the electron distribution. The quantity B c 
is then 


B c = 


_ 3.8 m e T l, e 


a e 


/l < g e,n v e 

V» (<r;v e > 


Rotational frequency . - From the imaginary part of equation (37), 

- i "R N <"i,e + 'V>- i f f “l +i f 5 - i |"l,l P=0 

an expression for u> R , the associated oscillatory frequency, can be obtained. 
o>j = 0, noting that /Lt ± . » /ij_ e > and simplifying give 

i * | [“K, i - ^1, i<L + M + P) - z] - P = 0 


which, when the Einstein diffusion relation is used, reduces to 


w R = - 
R B 


1 a?kT. 


Q_Q(L + M ).p 

LV e ^i,i/ N N 


From equations (12) and (13) the following expression can be obtained: 


■ + “V ‘ i. fm 0 = 0 

r 3r \ u 3r / 


D ' i ' B ' ~ 

r 3r \ 3r 


^ B criticaP 

(41) 

(42) 
Taking 

(43) 

(44) 

(45) 
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Integrating equation (45), setting the integration constant equal to zero, and noting that 
D l,i >>D l,e and ^1, i » ^1, e ^ve 



kT. 


_1 J. 

e n. 


0 



(46) 


Equation (46) can be used to evaluate the expressions for L, M, and P. Values of 

5 and ^ (L + M) + P appearing in equation (44) then become 
N N 


Q= 0.9 
N 


"\ 


> 


O 2 2kT i 

— (L + M) + P = 1.5 af. - 

N u e 


(47) 


Hence, from equations (47) and (40) and the fact that D, . » D, , equation (44) re- 

i e 

duces to the following expression for 


U R = 


kTj 

eB 


“ (°-9 a \ ~ 15 4) 


(48) 


Derivation of Ion Production Energy Per Beam Ion 

A classical derivation and a description of the anomalous aspects of the eV/ion for 
an ion thrustor are set forth in this section. The quantity eV/ion is defined as the power 
used in the plasma discharge divided by the beam current; hence, 


eV 

ion 


r e AV I 


(49) 


where I g is the cathode to anode current, AVj is the discharge voltage, and Ig is the 
ion beam current. The electron current density in the plasma source, with no instability 
assumed, can be written as 
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(50) 


J 


e 


= -D 


l, e 


e — + np, 
9r 


> e 


9V 

9r 


where the electrons are assumed to carry the current. Substituting the Einstein diffusion 
relation into equation (50) yields 


T _ ( kT e 1 9n 9V 

J e n/i l e " + — 

fc ,e \ e n 9r 9r . 


Using equations (16) and (46) and neglecting small terms give 


J e = n 0^i,e J l (o! 0 r ) eQ! 0 


(kT 


(51) 


(52) 


where J^ff^r) is a Bessel function of argument a^r. Finally, the current I g is equal 
to 2jralJ e (a), where a is the anode radius and SL is the length of the plasma. 

The beam current density Jg can be written as 

J B = env! (53) 

where v! is the velocity corresponding to one-half of the electron temperature. This 
assumption results from the Bohm criteria for stable plasma sheaths (ref. 11). The 

beam current is then equal to f 27rrJg dr. When equation (16) is used, the beam cur- 
rent becomes ® 


I B = en^vlca^ 


(54) 


where c is 2ir f - jQ(o'or)d( / -\. Taking the classical expression for 

** 0 a \a / ? 

p 

m e /eT e B gives the expression for the quantity eV/ion as 


as 


eV 

ion 


kT 

m e^l^ a 0 a ^ — “ 4. 8 nil AVj 
eT e v i c ®^ a ^ 


(55) 


where a^a is equal to 2. 4 as shown in equation (17). It can be seen that for this classi- 
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o 

cal derivation eV/ion behaves as l/B . 

In the case of anomalous diffusion, the current density in the plasma source can be 
represented by two terms (ref. 12 ): 


**e - ^class + ^anom 


(56) 


where J c j ass is the classical contribution and J 
The magnitude of with respect to J class 

netic field. Equation (56) can be rewritten as 


is the anomalous contribution. 

anom 

depends on the magnitude of the mag- 


J class= J e[ 1 -’< B >] 


(57) 


As the magnetic field becomes large, the quantity <J anom becomes much more significant 


than J 


class’ 


so that 1 - 77 (B) — 0 (ref. 12). The current density J can then be written 


from equations (57) and (51) as 



1 - 77(B) 


Following the foregoing derivation for eV/ion yields 


eV _ A 
ion ’ [l - T,(B)]B 2 


(58) 


(59) 


where A is a constant. From the condition that 1 - 77(B) — 0 for large magnetic field, 
an expansion of 1 - 77 (B) can properly be expressed as 


1 



(60) 


Thus, for the proper values of the constants A^, Ag, Ag, etc. , the value of eV/ion can 
vary as l/B (but not l/B 2 ), can be a constant, or can increase with increasing magnetic 
field. The results of using equation (56) and the condition 1 - 77 (B) — 0 for large mag- 
netic field are to show, mathematically, that the experimental results can be consistent 
with this explanation of anomalous diffusion. 
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RESULTS AND DISCUSSION 


A typical experimental variation of eV/ion with magnetic field is shown in sketch (b). 
The classical derivation given by equation (55) indicates that the variation of eV/ion with 

o 

magnetic field should behave as 1/B . Typical curves ap- 
pear to exhibit this behavior (ref. 13) at low magnetic field 
strengths. In the anomalous description of the eV/ion 
curve, it was shown that the curve could level off or in- 
crease at large values of magnetic field. This behavior is 
indicated by the right side of sketch (b) and by curves given 
in reference 14. The transition from the shape of eV/ion 
curves can be used to obtain experimental values of the 
critical magnetic field B . Values of the critical magnetic 

V 

field from equation (41) can thus be compared with these ex- 
perimental values. 



(b) 


Evaluation of Ion Production Energy Per Beam Ion at Critical Magnetic Field 

The magnitude of the quantity eV/ion at the critical magnetic field can be determined 
from equation (55) with B taken as B c determined from equation (41). (This quantity 
is evaluated in the next section. ) Selecting the following values of the variables of equa- 
tion (55) corresponding to a 10-centimeter-diameter ion thrustor with mercury as a pro- 
pellant 


a = — m B = 37 G (from eq. (41)) 

10 c 

Vj = 1. 4xl0 2 m/sec 
a = 0. 05 m 

ca 2 = 2vX5.4xi0“ 4 m 2 

yields a value of eV/ion equal to 320. Typical values of eV/ion for these conditions 
range from 400 to 600. The agreement is fair, which is all that should be expected from 
this derivation. 

The main objective of the derivation was to differentiate between the classical and 
anomalous behavior of the eV/ion against magnetic field curve. In achieving this, the be- 
havior of the eV/ion with respect to changes in magnetic field was of uppermost impor- 


AV = 50 V 


kT 


= 5 V 


T g = 6x10"® sec 
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tance. The derivation is not rigorous with respect to all the variables that might affect 

the eV/ion value of an ion thrustor. First, the blockage of the screens of the accelerator 

system is not taken into account since it is assumed in equation (53) that all the ions that 

reach downstream become part of the beam. Since the screens do block some of the ions 

from getting into the beam (physical blockage, about 50 percent with accelerating electric 

field effectively reducing this value), the calculated value of eV/ion is lower than the true 

value. Also, accelerating voltage and propellant utilization affect eV/ion. Neither of 

these factors is included in the foregoing derivation. It might be noted that thrustors 

having different anode radii have approximately the same value of eV/ion at B . The 

oo c 

classical expression, equation (55), predicts this fact since B*a is constant with re- 
spect to size, as can be seen from equation (41). This point will be given a more detailed 
explanation in the following section. 


Evaluation of Critical Magnetic Field and Oscillatory Frequency 


The equation for B c , equation (41), can be rewritten as 



(61) 


For a given propellant and with the electron distribution assumed not to change drasti- 
cally from one size ion thrustor to another, the right side of equation (61) is a constant. 
This condition implies that for the critical magnetic field B c the ratio of the ion thrustor 
radius to the Larmour radius is a constant. This condition is consistent with the scaling 
law derived from experimental performance of the electron-bombardment ion thrustor 
(ref. 3). 

Evaluating equation (41) for the critical magnetic field B c requires a knowledge of 
the electron distribution in an ion thrustor. Such information is quite scanty. One ap- 
proach that can be taken with respect to the electron energy distribution is to assume a 
Maxwellian distribution with an average energy of 5 electron volts, corresponding to the 
first excitation energy level of mercury. Equation (41) was evaluated for a 10-centimeter- 
diameter mercury ion thrustor with the result that B was 37 gauss. Typical values 
(refs. 13 and 14) of the critical magnetic field B c for the 10-centimeter-diameter mer- 
cury ion thrustor, evaluated from the transition region of eV/ion against magnetic field 
curves, range from 30 to 40 gauss. The agreement is thus quite good. As a further 
check, an experimental electron energy distribution for an ion plasma source with argon 
propellant, operating at the typical discharge voltage of 50 volts, was evaluated in refer- 
ence 6. Of the two distributions given in this reference, the one with an average energy 
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of 4 electron volts was selected to simulate a typical operating mercury bombardment ion 
thrustor operating at a discharge voltage of 50 volts. The value of velocity used in the 
expression for the Larmour radius corresponded to 4 electron volts, the average energy 
of the distribution. The value of B c for a 10-centimeter-diameter mercury ion thrustor 
was then calculated (by using this distribution) from equation (41), and it was found that 
B was 41 gauss. Finally, using the experimental electron distribution from reference 6 
corresponding to an average energy of 10 electron volts (as a maximum that might reason- 
ably be expected) and values of a g n and oj for argon yielded a value of B c from equa- 
tion (41) such that B was 23 gauss. From reference 12 the value of B„ for an argon 

c 

propellant ion thrustor, evaluated from the transition region of the eV/ion against mag- 
netic field curve, is between 20 and 30 gauss. 

To determine the critical magnetic field from equation (41), it is necessary to evalu- 
ate the behavior of the quantity \/(<t v ) /3 (cr!v ) . The fact that this quantity enters 
equation (41) as the square root diminishes the sensitivity of the critical magnetic field 
for changes of propellant as well as for corresponding electron distributions. Further- 
more, because the values of or! relative to a & for likely ion thrustor propellants do 
not differ greatly from one atom to another, the effect of propellant change on the critical 
magnetic field is not great. However, the electron distribution is a more important fac- 
tor in determining the quantity V( a e n v e^ 3 ^ a i v e^ ' T * ie electron energy distributions 
used in the calculations of B c had average energies of 4, 5, and 10 volts, corresponding 
to short to long tail types of electron energy distribution. The results were to obtain 
values of B c at 41, 37, and 23 gauss. It is reasonable to conclude that for a few typical 
propellants and for a variety of possible electron distributions the value of B for a 10- 
centimeter-diameter ion thrustor should fall between 23 and 41 gauss. This is consistent 
with the experimental values obtained from reference 14. 

The value of the oscillatory frequency associated with the anomalous diffusion can be 
calculated from equation (48). Considering a 10-centimeter -diameter engine with the 
critical magnetic field equal to 40 gauss and kT^/e equal to 1/20 electron volt, oj r is 
found to be 


co R = 2. 25X10 4 sec -1 

This is a frequency of 3500 cps. In reference 10 a plasma source with a magnetic field 
and a perpendicular electric field was described as operating in the region of ei hanced 
conduction. It was also found in this reference that a perturbation of number density 
existed, and associated with this perturbation was a frequency of approximately 
16 000 cps. Inasmuch as equation (48) is for the initiation of anomalous diffusion, rather 
than for an advanced stage of diffusion, a difference of less than an order of magnitude is 
felt to be substantial agreement. 
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CONCLUSIONS 


Expressions were derived for the onset of anomalous diffusion in an axially symmet- 
ric plasma with an axial magnetic field and a radially outward electric field. This 
plasma configuration corresponds to the ionization chamber of an electron-bombardment 
ion thrustor. A theory of plasma stability in a magnetic field first introduced by 
Kadomtsev and Nedospasov was modified and used in this derivation. As a result of this 
study, expressions for the critical magnetic field for the onset of anomalous diffusion and 
for the oscillatory frequency associated with the anomalous diffusion were obtained. The 
expression derived for the critical magnetic field varied inversely with anode radius in 
the same manner as an experimentally derived scaling law for the ion thrustor. Values 
of the critical magnetic field obtained from the derived expression compared well with 
experimental values obtained from curves of ion production energy per beam ion (eV/ion) 
against magnetic field. The oscillatory frequency associated with the anomalous diffusion 
fell within less than an order of magnitude of the frequency of oscillations found experi- 
mentally in a plasma source of similar configuration. Expressions for the variation of 
eV/ion with magnetic field were presented, covering both classical and anomalous re- 
gimes. In the classical region of this curve the value of eV/ion of the ion thrustor fell as 
l/B , where B is a magnetic field. It would be expected from these classical consider- 
ations that an increased magnetic field would contribute to exceptionally low values of 
eV/ion. It is postulated that this expectation is not realized because of the onset of anom- 
alous diffusion. The effect of this enhanced diffusion is to eliminate completely the effect 
of the magnetic field in lowering the eV/ion of the ion thrustor beyond the critical value 
of the magnetic field. In the anomalous region of the eV/ion against magnetic field curve, 
the variation is either constant or increasing with increasing magnetic field. This behav- 
ior corresponds to the increased diffusion of electrons in the anomalous diffusion regime. 
According to the applied theory, the onset of anomalous diffusion implies a perturbation 
of number density and potential with a 9 -direction dependence. As a result, a Hall effect 
causes increased conduction in the r-direction and, hence, an increase of current beyond 
the classical value calculated for the ion thrustor. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 30, 1966, 

120-26-02-10-22. 



APPENDIX A 


SYMBOLS 


A 

constant 

z 

z -direction in cylindrical coor 
dinates 

a 

anode radius 


B 

magnetic field 

“o 

2. 4/a 

C 

defined by eq. (34) 

a l 

3.8/a 

D 

diffusion coefficient 

77 (B) 

^anom/^E 

E 

electric field 

6 

9 -direction in cylindrical coor 

e 

electron charge 


dinates 

I 

current 

M 

mobility coefficient 

J 

current density 

V 

collision frequency 

J 0 

zero-order Bessel function 

<T 

collisional cross section 

J 1 

first-order Bessel function 


ionization cross section 

k 

Boltzmann constant 

T 

mean free collision time 

£ 

anode length 

0) 

oscillatory angular frequency 

m 

fundamental mode 

w e 

electron cyclotron frequency 

m e 

mass of electron 

Subscripts: 

n 

number density 

B 

beam 

r 

r-direction in cylindrical co- 

c 

critical 


ordinates 

e 

electron 

r T 

electron Larmour radius 



Jj 


I 

imaginary 

T 

temperature 

i 

ion 

t 

time 

n 

neutral 

V 

potential 

R 

real 

)— 1 
> 
<1 

discharge potential 

r 

radial direction 

V 

velocity 

z 

z -direction 

X 

argument of Bessel function 

e 

0-direction 

Z 

ion production coefficient 
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0 unperturbed 

1 perturbed 

1 perpendicular 


Superscript: 
— vector 


APPENDIX B 


EFFECTS OF CURRENTS IN z-DIRECTION 


If currents in the z-direction are taken into account, equations (4) and (7) become 


n v „ = -D 
e e, z € 


3n 


9z 


3n. 


e 3v 

+ — n c 


3z 




az * az 


(Bl) 

(B2) 


Since there is no applied electric field in the z-direction, and since it is assumed that the 
analysis is restricted to the central portion of the discharge far from both end walls, the 
potential gradient in equations (Bl) and (B2) can be neglected. As a result of equations 

p p 

(Bl) and (B2), equations (10) and (11) have the additional terms -D^ a n/3z and 
2 2 e 

-D^ 3 n/3z , respectively, on their right sides. As in reference 9, where a similar cir- 
cumstance is met, n(z) is assumed to vary as cos( 7 r/L z )z, where L z is some character- 
istic length and the unperturbed number density is then given as 

n Q = n(r)n(z) (B3) 

Equation (14) then has the additional term ^(j ir/L z )^ - D.ju^ e ( ff /L z )^JnQ on the 

right side. The constant A in equation (15) is now changed, but the solution of equa- 
tion (15) 


V r) = 4J 0 (<* 0 r) 


(B4) 


still remains intact. 

It is necessary to consider the effect of currents in the z-direction on the stability 
calculation. Equations (20) and (21) have the following additional terms on the right side: 



and 
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D 


e 




respectively. If, as before, the unperturbed potential gradient in the z -direction 
3Vq/3z is neglected for reasons given previously and if the perturbed quantities n^ and 
Vj are as given in this report, that is, as functions of r and 6 only, then all the addi- 
tional terms in this appendix drop out. Taking the perturbed quantities and as 
functions of only r and 6 is consistent with statements made in reference 9, where a 
similar situation is encountered. Thus, since the additional terms of equations (20) and 
(21) are dropped out, the analysis continues unchanged. 
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